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Abstract 

We characterize the biorthogonal polynomials that appear in the 
theory of coupled random matrices via a Riemann-Hilbert problem. 
Our Riemann-Hilbert problem is different from the ones that were 
proposed recently by Ercolani and McLaughlin, Kapaev, and Bertola 
et al. We believe that our formulation may be tractable to asymptotic 
analysis. 

1 Introduction 

The biorthogonal polynomials that appear in the theory of coupled random 
matrices E3 EU EE] are characterized by the property that 

J [ Pk{x)q {y)^ V{x) ' W{y)+2TXy dxdy = 0, if j + k, (1.1) 

where p^ and qj are polynomials of exact degrees k and j, respectively. In 
(jl.l|) we have that V, W : M — > M are given functions with sufficient increase 
at infinity so that the integrals converge, and r ^ is a nonzero coupling 
constant. The integration is over M 2 . 

Ercolani and McLaughlin [ IfJ showed that the two sequences of biorthog- 
onal polynomials (pk) and (qj) exist, that they are unique, and moreover, 
that pk has exactly k simple real zeros, see also |22j . They also gave a 
Riemann-Hilbert formulation for the biorthogonal polynomials which is non- 
local in character. Recently, for the case that V and W are polynomials, 
Kapaev JH] and Bertola et al. 0] gave local Riemann-Hilbert problems. If 
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d = deg W, then the Riemann-Hilbert problems for pk are formulated for 
d x d-matrix valued functions in |IJ ^] . 

In this note we derive a different Riemann-Hilbert problem. Our Riemann- 
Hilbert problem is based on the fact that the biorthogonal polynomials can 
be characterized as multiple orthogonal polynomials (see below). The for- 
mulation of a Riemann-Hilbert problem for multiple orthogonal polynomials 
is due to Van Assche et al. |21| 

An outstanding problem in random matrix theory is to provide a rigorous 
asymptotic of eigenvalue statistics for coupled random matrices. The basic 
example is the so-called 2-matrix model in which we have a probability 
measure on pairs (Mi, M 2 ) of Hermitian N x N matrices of the form 

— exp(-Tr(V(Mi) + W{M 2 ) - 2rM 1 M 2 ))dM 1 dM 2 . 
Zn 

Statistical quantities on eigenvalues of M\ and M 2 can be expressed in terms 
of the biorthogonal polynomials p^ and qj given by , see ^2 . The 
connection to biorthogonal polynomials would be very useful, if one has, in 
addition, a complete asymptotic description of the biorthogonal polynomi- 
als. Then it would be possible to compute eigenvalue statistices in the large 
N limit. Indeed, although the calculations are somewhat involved, this has 
been carried out in the Gaussian case V(x) = x 2 , W(y) = ay 2 , in 10:. Rig- 
orous asymptotics for biorthogonal polynomials with more general functions 
V and W are not known. 

In the 1-matrix case, the statistical quantities on eigenvalues are given 
in terms of orthogonal polynomials ^7], which have been characterized by 
a Riemann-Hilbert problem ^2]. The steepest descent / stationary phase 
method for Riemann-Hilbert problems was applied with great success to 
orthogonal polynomials [SlOilE]- As a result, the large N asymptotics of 
1-matrix models could be carried out in great detail, which in particular 
provided a proof of the universality of eigenvalue spacings for a large class 
of matrix models El Hfij - So there is hope that a similar asymptotic 
analysis of a Riemann-Hilbert problem for biorthogonal polynomials will 
lead to large N asymptotics for 2-matrix models. The formulation of a 
suitable Riemann-Hilbert problem is only a first step in this direction. 

For simplicity and clarity we formulate and prove our Riemann-Hilbert 
problem for the biorthogonal polynomial p^ for the first non-trivial case. 
This is the case where W is a polynomial of degree 4. Indeed, if W is a 
polynomial of degree 2, say W(y) = y 2 + 2by + c, then p^ is the orthogonal 
polynomial with respect to the weig ht e -V(x)+(rx-b) 2 Qn R and go there 

is a 

Riemann-Hilbert problem for p^ [7| I12j. 
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Thus we assume that W is a polynomial of degree 4 and we define 



y3 e -V(x)-W(y)+2rxy dy ^ ■ = ^ ^ % 



(1.2) 



These functions appear in the formulation of our Riemann-Hilbert problem. 
Riemann-Hilbert problem for Y 



The problem is to find a 4 x 4 matrix valued function Y : C \ R - 
having the following three properties. 

(a) Y is analytic on C \ 1. 

(b) Y has boundary values on R, denoted by Y + and YL, so that 



C 



4x4 



Y+(x) = Y-(x) 



(I w (x) wi(x) w 2 {x)\ 

1 



\0 



(c) Asz-* oo, we have 

Y(z) = ( / + 
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(1.3) 



(1.4) 



where k & No, n = [^] , n\ = [^] , and n 2 = [|] . Here [•] denotes 
the integer part. (Note that k = tiq + m + n.2-) 

The main result of this paper is that the Riemann-Hilbert problem for 
Y has a unique solution and that its (1,1) entry Y\\ is equal to the monic 
biorthogonal polynomial p/,- In what follows we use C(f) defined by 



C(f)(z) 



1 

2ttI 



■dx, 



zec\: 



x — z 

to denote the Cauchy transform of a function / : ]R — » R. 

Theorem 1.1. Lei f/ie functions Wj, j = 0,1,2, be given by (|1.2I) anc? 
/ef fc G Nq. Then the above Riemann-Hilbert problem for Y has a unique 
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solution given by 



Y 



( Pk 


C{pkW ) 




C{VkW2) \ 


(0) 

Pfc-l 




C(j>f_ x w x ) 




Pk-l 


C(P?2itoo) 






\pk-l 


C(p1-i^o) 







(1.5) 



where pk is the monic polynomial of degree k satisfying and £>jf2i? 

j = 0, 1, 2, are i/iree polynomials of degrees < k — 1. 

Remark 1.2. There is an immediate extension to polynomials of arbi- 
trary degree. If W is a polynomial of degree d, then the Riemann-Hilbert 
problem is for Y : C \ M -> C dxd so that 



y is analytic on C \ R. 



(1.6) 



The jump condition uses the d—1 functions uij(x) = J y^e v ( x ' w (y^ +2TX ydy, 
j = 0, 1, . . . , d — 2, and is given by 



y+(x) = Y-(x) 



(l w Q (x) w x (x) 

1 

1 

: : 




\0 

for x G R. The asymptotic condition is 



w d _3(a;) w d _ 2 (;c)\ 







1 / 



, (1.7) 



y(*) 



i + o 



diag z , z 



y /c no ni 



-1d-2 



(1.8 



where n, 



fc+rf-2-j 
<f-l 



for j = 0, 1, . . . , d — 2, and diag (•) denotes a diagonal 
matrix. 

The Riemann-Hilbert problem (|1.6|) - (|1.8|) has a unique solution and 
y.1 = Pk, where pk is the monic biorthogonal polynomial of degree k. The 
proof of the general case follows along the same lines as the proof of the case 
deg W = 4 that we will present in Sections 2 and 3 below. 
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Remark 1.3. By symmetry, there is a similar Riemann-Hilbert problem 
that characterizes the other biorthogonal polynomial qj, in the case that V 
is a polynomial. 

2 Multiple orthogonality 

We assume that W is a polynomial of degree 4. In this section we will 
characterize the monic biorthogonal polynomial pk of degree k through a set 
of orthogonality relations with respect to the three functions Wj, j = 0, 1, 2. 
As in Theorem 1.1, we use no = [^p] , ni = [^p] , and n-i = [§] . 

Lemma 2.1. We have 



Pk(x)x l Wj(x)dx = 0, for i = 0, 1, . . . , rij — 1, j = 0,1,2, (2.1) 



and these relations characterize the biorthogonal polynomial pk among all 
monic polynomials of degree k. 

Proof. Since W is a polynomial of degree 4, it is easy to see that 



where f^zi+j is a polynomial of exact degree 3i + j. For any function /, we 
then have if we integrate by parts i times 



If / is the biorthogonal polynomial p^ then the left-hand side of Q2.2JI is 
zero if 3i + j < k. This corresponds exactly with i < rij — 1 for j = 0,1,2, 
so that by the right-hand side we have the relations (|2.1|) . 

Conversely, if / is a monic polynomial of degree k that satisfies the 
relations J f{x)x l Wj(x)dx = for i = 0, . . . , rij — 1, j = 0,1,2, then the 





(2.2) 
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left-hand side of Q2.2|) is zero for % < rij — 1 and j = 0, 1, 2. The polynomials 
7T3j_|_j with i < rij — 1 and j = 0, 1, 2 are a basis for the polynomials of degree 
< k — 1. Hence / is the biorthogonal polynomial p k . □ 

Remark 2.2. The relations (|2.1j) are called multiple orthogonality relations 
of type II, see P EU [13 EH for more on this subject. 



3 Proof of Theorem 1.1 

Proof. We first establish uniqueness in the standard way. The only thing 
we have to observe is that both the jump matrix in (|1.3I) and the diagonal 
matrix in the right-hand side of (|1.4j) have determinant one. Then the proof 
of uniqueness follows as in [7[ Section 3.2]. 

We now prove that Y given by (|1.5|) satisfies the Riemann-Hilbert prob- 
lem. First we consider the first row of Y. The conditions (|1.3f) and (|1.4[) 
give for the (1,1) entry 

Y 11;+ = y n ,_, and Y u (z) = z k + 0(z k ^), 

These conditions are clearly satisfied if Y\i = pk, since pk is a monic poly- 
nomial of degree h. 

For the other entries in the first row, the jump condition (|1.3j) then is 

Y lji+ = Y^ _ + Y n -Wj-2 = Yij _ + PkWj-2 j = 2, 3, 4. 

By the Sokhotskii-Plemelj formula, this is satisfied by Yij = C{p k Wj-2)- 
The asymptotic condition (|1.4j) is 

Y lj (z) = 0(^- 2_1 ) asz^cx) (3.1) 

and we have to check that this is satisfied for Yij = C{p k Wj-2)- If we use 

1 "y^ 1 X i X n i~ 2 1 



X — Z ^ Z l+1 z n i~ 2 X — Z 

i=0 



then we see that for j = 2, 3, 4, 

1 f p k (x)Wj- 2 {x) 



C(p k Wj- 2 )(z) 



2m J x 



^2 j Pk(x)x l Wj-2(x)dx^ z % 1 

1 f p k (x)wj. 2 (x)x n ^ \ 2 



+ [— r^">-3-^j- dx 2 -n,_ 2 _ (32) 

1 2m / x — z 
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Because of the multiple orthogonal relations (|2.1j) satisfied by pk, ()3.2[) re- 
duces to 

ccp^w = (- y ^ — a*) 

which shows that (|3.1|) is indeed satisfied if Yy = C{pkWj-2) 

Next we consider the second row of Y . The conditions Y21 + = ^21,- an d 
Y2i(#) = 0(z fc_1 ) as 2 — > 00 are clearly satisfied if Y21 is a polynomial 
of degree < k — 1. The jump conditions for the other entries in the second 
row 

Y 2j ,+ = Y 2j - + = Y 2 j- +pf_ 1 w j - 2 

are then also satisfied if Y 2 j = Cip^^Wj^) for j = 2,3,4. We need to be 

able to choose p^_i so that the asymptotic condition (|1.3|) is also satisfied, 
which means that 

C(p { u\w ) = z~ n * +O(z- no - 1 ), 

(3.3) 

C(Pfc-i«>i) = O^-^'" 1 ) for j = 1, 2. 

Expanding C(p^ 1 Wj) as in (|3.2j) . we see that (|3.3(l is satisfied if p^f^ is such 
that 

pfl 1 (x)x l wo{x)dx = 0, for i = 0, . . . , no — 2, (3.4) 
p^ 1 (x)x' l wo(x)dx = —2iri, for i = no — 1, (3-5) 
p^ 1 (x)x i tt>j(x)(ia; : 

The conditions Q3.4JI and (|3.6j) give fc — 1 homogeneous conditions on the 
k free coefficients of and so there exists a non-zero polynomial 

satisfying these conditions. To be able to have Q3.5|) as well we must exclude 
the possibility that 



0. 


for i 


= 0,. 


— 2ni, 


for i 


= n 


0. 


for i 


= 0,. 



pf_ l {x)x n °- l w {x)dx = 0. (3.7) 

However, if Q3.7[) would hold, then pp. + would be a monic polynomial 
of degree k that satisfies the multiple orthogonality relations (|2.1j) . which 
is impossible, since these relations characterize the biorthogonal polynomial 
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Pk by Lemma 2.1. Thus ()3.7j) cannot hold. Then we can normalize by 
multiplying it with a suitable constant, so that q3.5|) is satisfied. This proves 

that we can indeed choose pi_i so that the second row of V satisfies all the 
conditions imposed by the Riemann-Hilbert problem. 

In exactly the same way, we handle the third and fourth rows of Y. 

This completes the proof of Theorem 1.1. □ 

4 Conclusion 

We have characterized the biorthogonal polynomials that appear in the the- 
ory of coupled random matrices via a Riemann-Hilbert problem which is 
differnt from the Riemann-Hilbert problems derived in [1J ^3 . Recent expe- 
rience [13 with similar higher order Riemann-Hilbert problems leads us to 
believe that our Riemann-Hilbert problem may be tractable to asymptotic 
analysis. However, up to now we have not been able to apply the steep- 
est descent method succesfully to this problem, and the actual asymptotic 
analysis of biorthogonal polynomials remains a major open problem. 

References 

[1] A. I. Aptekarev, Multiple orthogonal polynomials, J. Comput. Appl. 
Math. 99 (1998), 423-447. 

[2] A. I. Aptekarev, A. Branquinho, and W. Van Assche, Multiple orthog- 
onal polynomials for classical weights, Trans. Amer. Math. Soc. 355 
(2003), 3887-3914. 

[3] M. Bertola, B. Eynard, and J. Harnad, Duality, biorthogonal polynomi- 
als and multi-matrix models, Comm. Math. Phys. 229 (2002) 73-120. 

[4] M. Bertola, B. Eynard, and J. Harnad, Differential systems for 
biorthogonal polynomials appearing in 2-matrix models and the as- 
sociated Riemann-Hilbert problem, preprint nlin.SI/0208002 

[5] P. Bleher and A. Its, Semiclassical asymptotics of orthogonal polyno- 
mials, Riemann-Hilbert problem, and universality in the matrix model, 
Ann. Math. 150 (1999), 185-266. 

[6] P. M. Bleher and A. B. J. Kuijlaars, Random matrices with exter- 
nal source and multiple orthogonal polynomials, to appear in Internat. 
Math. Research Notices (2004), preprint math-ph/0307055 . 



8 



[7] P. Deift, Orthogonal Polynomials and Random Matrices: a Riemann- 
Hilbert Approach, Courant Lecture Notes in Mathematics, 3. New York 
University, New York, 1999. Reprinted by Amer. Math. Soc, Provi- 
dence R.I., 2000. 

[8] P. Deift, T. Kriecherbauer, K. T-R McLaughlin, S. Venakides, and X. 
Zhou, Uniform asymptotics for polynomials orthogonal with respect to 
varying exponential weights and applications to universality questions 
in random matrix theory, Comm. Pure Appl. Math. 52 (1999), 1335- 
1425. 

[9] P. Deift, T. Kriecherbauer, K. T-R McLaughlin, S. Venakides, and X. 
Zhou, Strong asymptotics of orthogonal polynomials with respect to 
exponential weights, Comm. Pure Appl. Math. 52 (1999), 1491-1552. 

[10] N. Ercolani and K. T-R McLaughlin, Asymptotics and integrable struc- 
tures for biorthogonal polynomials associated to a random two-matrix 
model, Physica D 152/153 (2001), 232-268. 

[11] B. Eynard and M. L. Mehta, Matrices coupled in a chain: eigenvalue 
correlations, J. Phys. A 31 (1998), 4449-4456. 

[12] A. S. Fokas, A. R. Its and A. V. Kitaev, The isomonodromy approach to 
matrix models in 2D quantum gravity, Comm. Math. Phys. 147 (1992), 
395-430. 

[13] A. A. Kapaev, The Riemann-Hilbert problem for the bi-orthogonal 
polynomials, J. Phys. A 36 (2003), 4629-4640. 

[14] A. B. J. Kuijlaars, K. T-R McLaughlin, W. Van Assche, and M. 
Vanlessen, The Riemann-Hilbert approach to strong asymptotics for 
orthogonal polynomials on [-1,1], to appear in Adv. Math, preprint 
math.CA/0111252[ 



[15] A. B. J. Kuijlaars, W. Van Assche, and F. Wielonsky, Quadratic 
Hermite-Pade approximation to the exponential function, a Riemann- 



Hilbert approach, preprint math.CA/0302357 



[16] A. B. J. Kuijlaars and M. Vanlessen, Universality for eigenvalue cor- 
relations from the modified Jacobi unitary ensemble, Internat. Math. 
Research Notices 2002 (2002), 1575-1600. 

[17] M. L. Mehta, Random Matrices, second edition, Academic Press, 
Boston, 1991. 



9 



[18] M. L. Mehta, A method of integration over matrix variables, Comm. 
Math. Phys. 79 (1981), 327-340. 

[19] E. Nikishin and V. Sorokin, Rational Approximation and Orthogonal- 
ity, Translations of Mathematical Monographs 92, Amer. Math. Soc. 
Providence R.I., 1991. 

[20] W. Van Assche and E. Coussement, Some classical multiple orthogonal 
polynomials, J. Comput. Appl. Math. 127 (2001), 317-347. 

[21] W. Van Assche, J. Geronimo and A. B. J. Kuijlaars, Riemann-Hilbert 
problems for multiple orthogonal polynomials, in: NATO ASI Special 
Functions 2000. Current Perspective and Future Directions (J. Bustoz 
et al. eds.), Kluwer Academic Publishers, Dordrecht, 2001, pp. 23-59. 

[22] Y. Xu, Biorthogonal polynomials and total positive functions, J. Phys. 
A 35 (2002), 5499-5510. 

A. B. J. Kuijlaars 

Department of Mathematics 

Katholieke Universiteit Leuven 

Celestijnenlaan 200B 

B-3001 Leuven, BELGIUM, 

E-mail: arno@wis.kuleuven.ac.be 

K. T-R McLaughlin 
Department of Mathematics 

University of North Carolina at Chapel Hill 
Chapel Hill NC 27599, U.S.A. 
E-mail: mcl@amath.unc.edu 



10 



